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O ■ ABSTRACT 

\o ■ 

, Context. A specific example of Bianchi Type VII;, models, i.e. those including universal rotation (vorticity) and differential expansion (shear), 
• has been shown in Jatfe et al. (2005) to correlate unexpectedly with the WMAP first-year data. 
' Aims. We re-assess the signature of this model in the WMAP 3-year data. 
Methods. The cross-correlation methods are described in Jalfe et al. (2006a). We use the WMAP 3-year data release, including maps for 
O I individual years, and perform additional comparisons to assess the influence of both noise and residual foregrounds and eliminate potential 
non-cosmological sources for the correlation. 

Results. We confirm that the signal is detected in both the combined 3-year data and the individual yearly sky maps at a level consistent with 
^ our original analysis. The significance of the correlation is not alfected by either noise or foreground residuals. 
Conclusions. The results of our previous study are unchanged. 
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. . . 1 . Introduction pipelines, calibration and corrections for foreground contami- 
nation. 

In Llaffe et alJ ll2005[l2006all . we reported on the unexpected de- 
tection of a correlation between the cosmic microwave back- 2. Methods 

ground (CMB) sky measured by the WMAP first-year data re- ■ ■ ■ . 

lease and a Bianchi Type VII,, template. This result is partic- The methods used are described in detail in| Jaffe et alj (|2006a|). 

ularly provocative in that such a signal may provide an ex- I" brief, we use a minimum-^^ j^gj^od to estimate the cross- 

planation for several indep endendy discovered ano malies: the correlation amplitude between the WMAP data and thedeter- 

low quadrupole amplitude llde Oliveira-Costa et alJl20 04) rela- mini^tic Bianchi template derived by BarrowetalJ USM and 

live to the best-fit cosmological model, the curious alignment extended in| Jaffe et al. ,a006b) to include arbitrary dark energy 

of several larg e angular scale multipole s along the so-called ^nodels. 

"axis of evil" dl^and & MagueiiolEooi . a sign ificant power F^^H ^ky analysis is only possible using the highly pro- 
asymmetry between two hemispheres on the sky ("H ansen et al.l messed internal linear combination (ILC) map, created by com- 
[2004; Eriksen et al. 2004b) and a large temperature decre- fining the data at diff^erent frequencies using weights that min- 
ment towards a specific part of th e south ern sky as re vealed imize the variance while preserving the cosmic signal. Such a 
by a wavelet analysis (ICruz etalJ Eooi IVielva et al.' "2004). combination is expected to clean the data of galactic emission. 
The subtraction of the Bianchi component leaves a statistically However, as pointed out by [ Eriksen et al. |(|2004a), the method 
isotropic sky does leave significant foregro und residuals which particularly 

affect the galactic plane region. Einshaw etal. (.2006.) point out 

In this n ote, we examine the ne wly released three-year data that this is a consequence of the method minimising the vari- 

described in ' Hinshaw et all (''2006) and show that the signal is ance by anti-correlating the CMB and foregrounds. In the 3- 

detected both in the combined 3-year data and in each indi- year release, an additional "de-biasing" is therefore applied that 

vidual yearly map, despite differences in the data processing attempts to correct for this. With these data, we use the ef- 
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ficient total convolver algorithm jWandelt & G6rskill200ll) in 
harmonic space to find the best orientation of the Bianchi tem- 
plate relative to the data. We then use the sky maps from ten in- 
dividual differencing assemblies (DAs) spanning five frequen- 
cies (K, Ka, Q, V, W) to confirm this result and assess the extent 
to which the full sky fit is affected by foreground residuals. 

For individual frequency bands (derived by co-adding the 
corresponding DAs together) and combinations of bands, we 
fix the template at its best-fit location and perform simultane- 
ous fits in pixel space of the Bianchi template and three tem- 
plates to model the galactic foreground contribution. For ther- 
mal dust emission, we use either the lOOjum sky map derived 
by Schlegel et al. (1998, hereafter SFD) or a template at an 
eff'ective frequency of 94 GHz (the WMAP W-band) derived 
from SFD using model 8 of Finkbeiner et al. (1999, hereafter 
FDS) that includes an additional correction for the dust tem- 
perature. Free-free emission is traced by maps of H„. We con- 
sider two templates that differ in the way that they combine 
the available H p, survey da t a - D ickinson et al. (2003, here- 
after DDD) or iFinkbeine? feOOS). For Galactic synchrotron 
emission we a gain consider two alternative templates: the 408 
MHz survey of lHaslam et al.l(ll982h . and, following the WMAP 
analysis, the K-Ka diff'erence map. The advantage of the lat- 
ter is that it will more realistically account for the synchtron 
emission at the frequencies of interest, capturing some vari- 
ation in the synchrotron spectral behaviour that a simple ex- 
trapolation of the 408 MHz data cannot. Moreover, it will also 
trace any additional foreground component at low frequencies 
that is not represented by the dust or H„ templates. This in- 
cludes at least some fraction of the co-called "anomalous dust" 
emission - due to its correlation with the dust template - al- 
though the physical origin of this component still remains un- 
clear (see Davies et al. 2006 for a relevant discussion). The 
fits are confined to the region outside the conservative KpO 
mask that eliminates bright emission from the Galactic plane 
and several hundred point sources. This is necessary since the 
former exhibits complex spatially varying spectral behaviour 
that the templates are not sufficiently accurate to quantify. We 
estimate the effect of high-latitude foreground residuals not 
well traced by these templates by fitting the Bianchi mo dels 
to the foreground maps derived in i Hinshaw et all ( l2006l) us- 
ing a maximum-entropy (MEM) method. These include pixel- 
by-pixel variations in the foreground emission, but their noise 
properties are usual considered too complex to include in a cos- 
mological analysis. However, this is not an issue for this work. 

In order to characterize the effect of noise and systematic 
artifacts that are not well described by the WMAP noise model, 
we utiUze year-to-year difference maps. 

Finally, we use the LILC simulation pipeline of 
lEriksen et alJ ( l2004al) to characterize the null hypothesis, 
namely the distribution of best-fit chance alignments arising 
between the Bianchi template and simulated skies containing 
no such contribution. A significance measure for a given tem- 
plate is then determined from the number of chance alignments 
with higher amplitude than that measured using the data. 
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Fig. 1. Significance contours as percentage of LILC simula- 
tions whose best-fit chance alignment amplitude is lower. The 
left (right) column shows the left-handed (right-handed) mod- 
els. Over plotted contours are at 99.3 and 99.5%. Two color 
scales are used to show the global structure (top) as well as that 
near the peaks (bottom). 

For computational efficiency, all maps are smoothed to a 
5°. 5 FWHM beam and downgraded to HEALPix' N^ide = 32 
resolution. All pixel-space fits include simultaneous fits of the 
monopole and dipole components, since these may contain ad- 
ditional non-cosmological contributions. 

3. Results 

The results of the model-space search usin g the ILC map ove r 
the full sky are presented in Fig. [2 As in ljaffe et al.l (l2006al) . 
there is a region of high significance in the parameter space 
for right-handed models around {x, Qo) = (0.6, 0.5) seen in the 
right-hand column of Fig. In the previous work, there was 
also a possible detection of a left-handed model at (x, Qq) - 
(0.6, 0.15), which we examined in detail as well and concluded 
was due to foreground contamination in the galactic plane re- 
gion. The left-hand column of Fig. [2 shows that this model is 
not a significant fit to the 3-year data. This change is likely due 
to the differences introduced by the de-biasing applied to the 
ILC in the 3-year release. This analysis repeats the search for 
the best-fit orientation of the template relative to the data and 
finds the same position as the previous best-fit to within one bin 
of 2.8° size at this resolution. Using the full-sky ILC maps, the 
best-fit model amplitude drops a few percent from 4.33 x 10"'" 
to 4.19 X 10"'", which corresponds to a change in significance 
from 99.8% to 99.5%. 

Fig-Elpresents several difference maps that indicate the am- 
plitude and morphology of noise, calibration, and map-making 
effects which might influence the analysis. Fig. |2 (a) shows 

' http://healpix.jpl.nasa.gov/ 
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Fig. 2. Difference maps showing noise and calibration effects. 
The monopole and dipole terms are removed. See text for de- 
tails. 



the Ka-band difference between the 3-year and 1-year map, 
where the time variability in a few point sources (shown here 
highly smoothed) is the dominant component. Fig.|2(b) depicts 
the differences in the P' year Ka-band sky maps arising from 
changes in the map-making and calibration methods between 
the current and original data releases. The observed structure 
arises from changes in the gain calibration solutions, and is 
particularly notable in the smoothed diffe rences ar ound point- 
sources along the galactic plane. See Jarosik et alJ ( l2006i) for a 
discussion of these differences. Figs.|2](c) and (d) are the same 
for the W-band, although the point sources are no longer a sig- 
nificant contaminant. Figs.|2(e) and (f) show two examples of 
the noise differences in the W-band between pairs of individual 
years. Fig.|2(g) gives the residuals from fitting the foreground 
templates to the MEM Q-band foregrounds (see below). Lastly, 
Fig.|2|(h) shows the Bianchi template at one quarter of its best- 
fit amplitude for comparison. The Bianchi anisotropy ampli- 
tude is clearly much higher than any of these effects. 

In Fig. |3l we plot the fits to year-to-year difference maps 
for each DA. The scatter clearly increases at high frequency, 
as expected due to the increased noise. W4 in particular shows 
significantly higher scatter than the other W-band assemblies, 
which is unsurpr ising considering th at it has the highest level 
of 1// noise (see Llarosik et al.ir2003l) . We have also performed 
fits to maps constructed from the MEM foreground solutions to 
test for the effect of residuals not traced by the foreground tem- 
plates. We find that the maximum Bianchi amplitude as fitted to 
these maps is less than 3% of the ILC fit amplitude, highest for 
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Fig. 3. Fits to year-to-year difference maps by assembly. The 
different symbols show different pairs of years, while the color 
indicates the differencing assembly. The filled black symbols 
are the averages over the different assemblies. 



the Q-band fits as expected for foregrounds that increasingly 
dominate at lower frequencies. The fits are plotted in Fig. |3 
and the grey bands reflect the uncertainty in the fit results due 
to the noise and foreground residuals. The error bars are the sta- 
tistical uncertainties almost entirely due to the isotropic CMB 
contribution. 

Tabled shows that the results of the Bianchi fitting to the 
WMAP data are very stable and largely robust to foreground 
correction method and template selection. Although the fits 
to individual frequency bands are slightly lower than the ILC 
fit amplitude, they remain near or above the 99% significance 
level. The most consistent results are derived when the K-Ka 
map is used as a synchrotron tracer, presumably because it can 
also account for whatever low-frequency foreground emission 
is not traced by the dust and free-free templates. By compari- 
son, using the template combination of EDS dust, Finkbeiner 
Ho, and Haslam for synchrotron (as in our original analy- 
sis), there remains some low-frequency residual which anti- 
correlates with the CMB and therefore lowers the fit ampli- 
tudes. Using the DDD instead of the Finkbeiner template 
and the SED instead of the EDS thermal dust template again 
yields very consistent results. There is a small difference at W- 
band where the SED template perhaps does not account for the 
thermal dust emission as well as the EDS template does. 

The difference maps, which have no cosmic signal, give fit 
amplitudes of only 3% or less of the signal, consistent with the 
estimate given by the MEM fits. The 3-year fit values are also 
plotted in Fig. |4] with fits to noise and foreground maps ad- 
ditionally shown for comparison. The difference between the 
3-year fits and the individual years is too small to see on this 
plot. The grey band around the fit value obtained with the ILC 
shows that, though they are slightly lower, the cut-sky fits do 
not vary in excess of what is expected due to such small fore- 
ground residuals and noise. Even with these effects, the fits are 
consistently at the 99% significance level compared with sim- 
ulations (shown by the dashed line). 
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Table 1. Amplitudes (i.e. the shear value (cr/H)n x 1 ) of the best-fit model derived from various combinations of data and 
various methods as described in the text of Jaff e et al.lll20 06a'). Uncertainties are shown for the ILC only, and they are roughly the 
same for all, since the noise is negligible for this low-resolution analysis and the structure outside the Galactic cut. For the WILC 
maps, the total convolver method is used on the fuU sky. For the individual bands, the KpO mask was imposed and foreground 
templates fit simultaneously for the remaining maps: in columns labeled ( i), we used the FDS dust, Finkbeiner Ha, and Haslam for 
synchrotron; in columns labeled ( ii), we used WMAP K-Ka for synchrotron; in columns labeled ( Hi), we use alternate templates 
SFD for dust and DDD for H^, with again WMAP K-Ka for synchrotron. The last columns show the significance measure for the 
3-year results, i.e. the percentage of simulations with lower amplitude. The significances of the difference map fits are not quoted 
since they would reflect the fact that there are more residuals in the data than in idealised simulations where the fitted templates 
are themselves used to generate the foreground contribution. 
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Fig. 4. WMAP 3-year fit amplitudes for the three different sets 
of templates labeled columns /, //, and //; in Tabled Also plot- 
ted are fit results from difference maps between years and from 
the MEM foreground maps. These indicate the degree to which 
noise and foreground residuals may be contaminating the fits, 
shown as the grey band. This band is also shown around the 
best-fit ILC amplitude for comparison. The dashed and dot- 
dashed lines represent the approximate amplitudes correspond- 
ing to 99% and 95% significance, respectively. 

4. Conclusions 

In I.Taff^eetal.1 ll200.5[ l20n6ah . we showed that removing a 
Bianchi component from the WMAP initial data release can 
account for several large angular scale anomalies and yield a 
corrected sky that is statistically isotropic. Although we do not 



repeat all of these analyses here, we simply note that the change 
in the fit amplitude is too small to alter the previous conclu- 
sions. We have additionally shown in this work that the correla- 
tion cannot be attributed either to instrument noise, foreground 
residuals, or any known systematic corrected in the three-year 
WMAP data release. 

In Jaffe et al. ( 2006b), we found that the cosmological pa- 
rameters required to reproduce the Bianchi morphology are not 
compatible with the "cosmic concordance" of other datasets, 
including the small-scale CMB structure itself, ruling out these 
models as a realistic physical explanation for the anomalies. 
Nevertheless, a recent analysis with a statistically rigorous 
method in Bridaes et al. ( 2006) has demonstrated, in agreement 
with our results, that the data do appear to require the Bianchi 
template. They give the "Bayesian evidence" in favor of a cos- 
mology including vorticity and shear compared to one without 
as roughly AlnE ~ 2, which is considered "significant". 

We conclude, therefore, that in the absence of an unknown 
systematic effect which could explain both the anomalies and 
the correlation, the WMAP data require an addition to the stan- 
dard cosmological model that resembles the Bianchi morphol- 
ogy- 
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